Order parameter in superconductors with non-degenerate bands 
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In noncentrosymmetric metals, the spin degeneracy of the electronic bands is lifted by spin-orbit 
coupling. We consider general symmetry properties of the pairing function A(k) in noncentrosym- 
metric superconductors with spin-orbit coupling (NSC), including CePtaSi, Ulr and Cd2Re2 07. We 
find that A(k) — x(k)t(k), where x(k) is an even function which transforms according to the ir- 
reducible representations of the crystallographic point group and t(k) is a model dependent phase 
factor. We consider tunneling between a NSC and a conventional superconductor. It is found that, 
in terms of thermodynamical properties as well as the Josephson effect, the state of NSC resembles 
a singlet superconducting state with gap function x(k). 

PACS numbers: 74.20.-z, 74.20.Rp, 74.50.-fr, 71.27.-fa 



I. INTRODUCTION 

A better understanding of unconventional supercon- 
ductivity has been sought ever since the pioneering dis- 
covery of heavy fermion superconductivity in CeCu2Si2ii 
Soon after, the effect of spin-orbit coupling (SOC) was 
recognized to play an important role in the superconduct- 
ing (SC) properties of heavy fermion materials, because it 
breaks spin rotation symmetry^i^ If there is no spin-orbit 
coupling, or if the crystal has space inversion and time re- 
versal symmetry, the electronic bands are doubly degen- 
erate everywhere in the Brillouin zone.- Then, the super- 
conducting states can be classified as spin-singlet (even 
parity) and spin-triplet (odd parity), where the spin is 
replaced by pscudospin when SOC is present iS^^iSii How- 
ever, if inversion symmetry is absent and SOC is present, 
the degeneracy of electronic band states is removed at all 
points in the Brillouin zone except at some highly sym- 
metric positions, while time reversal symmetry ensures 
that states with opposite momenta are degenerated 

The recent interest in noncentrosymmetric supercon- 
ductors was stimulated by the discovery of superconduc- 
tivity in CePtsSi and UlriSiS SOC is normally strong 
in cerium and uranium compounds since the atoms are 
heavy. A recent band structure calculation of CePtaSi 
found that the bands, which would be degenerate if SOC 
was absent, are split by 50-200 meV for the states close 
to Fermi leveliiS This splitting energy is a factor of more 
than a thousand on the characteristic energy scale ksTc. 
Therefore, the one-band theory, which applies to the 
spin-degenerate case, should be reformulated^ 

Actually, several superconducting materials in which 
inversion symmetry is absent have been known for some 
time. Sesquicarbide materials, with chemical formula 
R2C3_y (R = La or Y, which can be partially substi- 
tuted by a number of elements) have Tc's of up to 18 
K, and space group symmetry /43dii2iiii4 This belongs 
to the tetrahedral crystallographic class Td which has no 
symmetry centre. However, in sesquicarbides SOC does 
not seem to be important for conduction electrons, i.e. 
the bands are still spin degenerate >i^ii& 



In 2001, superconductivity with Tc = 1 K was re- 
ported in Cd2Re207iiii^ At room temperature, this 
material has the ideal pyrochlore structure, which in- 
cludes a centre of symmetry. However, superconductiv- 
ity occurs on the background of a noncentrosymmetric 
tetragonal crystal field after a series of structural phase 
transitions>i2iSSi2i In addition, SOC dramatically affects 
the electronic band structure, as shown by first-principles 
calculations for the pyrochlore phase i^^i^*^ Following fur- 
ther calculations)^ it was predicted that if the centre of 
symmetry is removed from the structure, the spin-orbit 
splitting of the bands reaches 68 meV ~ 700 ksTc. This 
result was indirectly verified by a photoemission study, 
which found that the energy spectrum is shifted notice- 
ably towards higher binding energies at low tempera- 
tures as compared to room temperature^ Therefore, 
Cd2Re207 should be considered in the same category as 
CePtsSi and Ulr. Moreover, Cd2Re207 is a useful test 
system for the theory since its superconducting proper- 
ties are well established JS^SS; The theory must, therefore, 
include the possibility for a nodeless s-wave-like order 
parameter!^ Another appealing feature of Cd2Re2 07 is 
the relative simplicity of its normal metallic state: it nei- 
ther orders magnetically)2& nor demonstrates Kondo-like 
behaviour as that found in CePtaSii^ 



In this Article, we consider general symmetry proper- 
ties of the superconducting state in noncentrosymmet- 
ric crystals in which the degeneracy is lifted by strong 
SOC. As a first approximation, we restrict consideration 
to a one-band model, neglecting possible pairing in other 
bands and inter-band interactions. Sec. ^discusses the 
symmetry properties of the superconducting order pa- 
rameter. In Sec. mil we consider Josephson tunneling 
between a conventional superconductor and a supercon- 
ductor with nondegenerate bands. Sec. IIVI discusses the 
limiting case of small SOC, using the Rashba Hamilto- 
nian as a specific example. In Sec. we discuss the 
results of our approach in relation to previous theoreti- 
cal developments and the applicability to real materials. 



II. SUPERCONDUCTING ORDER 
PARAMETER 

Throughout this paper we use the weak couphng ap- 
proach to treat the pairing interaction. While this ap- 
proach offers no insight about the pairing mechanism, it 
has the advantage that many resuhs can be obtained ex- 
actly. Also, the symmetry properties of the SC state re- 
main essentially the same as in strong coupling models.^ 



A. Single-particle Hamiltonian 

In the condensed state it is normally sufficient to treat 
relativistic effects by the Pauli approximation to Dirac's 
fully relativistic approach. We begin by briefly reviewing 
the properties of the single-particle Hamiltonian^ 
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where p = — iV is the momentum operator in the coordi- 
nate representation, ^(r) is the periodic potential of the 
crystal lattice, and <t = {crx,(Ty,(Tz) are the Pauli matri- 
ces. Here and in the following we use the unit system 
in which h = ks — ^- The space group of the crystal is 
defined as the set of operations g (rotations, reflections 
and translations) acting on the real space coordinates r 
which leave ^(r) invariant. 



g : V(r) = V{gr) = T/(r). 



(2) 



Because of the last term in (QJ, the transformation prop- 
erties of Hi are expressed as 



g : ffi(r) = V{g)Hi{gr)V{g)-' = Hi{r), 



(3) 



where ^{g) is a 2 x 2 spin-1/2 rotation matrix, ^{g) is 
defined as follows. If 5 is a bare rotation about the axis 
n by an angle </>, then 



T^ig) = cos in- (T sin — . 



(4) 



If g includes inversion I, g — Ig' , or a translation t(R), 
g = T(R)g', where g' is a bare rotation, then ©(g) = 
I? (5'), since inversion and translations do not change the 
spinor components. 

The solutions $(r) of the Shrodinger equation 

Fi*(r) = ^^(r) (5) 

are two-component spinors which transform according to 



g : *(r) = V{g)^igr) 



(6) 



In nonmagnetic crystals V{r) does not depend on spin, 
and the Hamiltonian Q is also invariant under time re- 
versal /C, 

/C : i/i(r) = {-iay)H*i{r){^i<jy)-' = i7i(r), (7) 



which uses the fact that spin reverses sign under time 
reversal, expressed by ayiJ*<Ty — —or. Correspondingly, 

K. : ^{v) = i-iayWiv). (8) 

Further, if t(R) is a proper lattice translation, 

T(R):i?i(r)=i/i(r + R)=i7i(r). (9) 

Therefore, as well as in the case of vanishing SOC, the 
solutions of the Shrodinger equation (|SJ| are Bloch waves 
labeled by a quantum number quasimomentum k, 



*k(r) = [/k(r)e'''^ 
where f7k(r) is the periodic spinor 



(10) 
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[/k(r + R) = t/k(r). (11) 



It follows from Eq. © that 

5:vl/k(r) = I?(g)[/k(gr)e^('''«"^) 
= V{g)U^{gv)e^^^'''^' 



■)^*(^_\,(r). (12) 



,(9) 



Clearly, ^g-ikli") is a solution of ©, corresponding to 
the quasimomentum g~^k. On the other hand, the same 
can be said about the function ^g-ik(r) obtained from 
^k(r) by the bare replacement k -^ 5^^k. If the spin 
degeneracy of the bands is lifted, these two functions can 
only differ by a phase factor that does not depend on r. 



■^[td^) = exp[ja3(k)]vE'^-,k(r). 



(13) 



This argument also applies to time reversal operation. 



/C:vl/k(r) ^ {-zay)U^{r)e 



ikr _ iTfCC) 



*L\^(r). (14) 



Let us define i(k) as the phase factor in 

vI/L^j,)(r) = i(k)vl/_k(r), t*{k) ^ t"i(k). (15) 

Using the antilinear property of the time reversal opera- 
tor /C : a^k(r) = ci*IC : ^k(r) and the fermionic nature 
of ^k(r), one obtains, 



/C^ : ^k(r) 



-*k(r) 

t*(k) /C : *_k(r) 

t*(k)t(-k)*k(r) 



Thus, i(k) is an odd function of k, 
t{-k) = -i(k) 



(16) 



(17) 



Before discussing the interaction part of the Hamilto- 
nian, we formulate the symmetry properties of the Hamil- 
tonian (^ using t he language of second quantization [See, 
e.g., Refs. !28l23|. The A^-particle wave function ^n is 
the anti-symmetrized product of N single-particle func- 
tions ^k(r) taken at N different points k. By writing 



(. . . , Ik, . . .) as the argument of ^n we emphasize that 
it describes a state in which a single-particle state with 
quantum number k is filled. By definition, the creation 
operator cj^, 
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.)=±$Ar(...,lk,...), 
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where the sign is defined by the antisymmetrization. In 
particular, we note that the matrix elements of the cre- 
ation operator are real numbers. Using the definition H18|) 
and the transformation properties of ^k(r), we obtain^ 
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(19) 



where Ui{9) is a gauge transformation. The transforma- 
tion properties of Ck follow from the fact that the number 
operator cj^Ck must be a scalar. 

The one-band single-particle Hamiltonian can be ex- 
pressed as 



Wi = X! ^k4ck. 



(20) 



where £,k — e^ — ^, Sk is the eigenvalue of the Hamil- 
tonian 1^ corresponding to momentum k, and /i is the 
chemical potential. We note that c^ and Ck are both two 
component spinors, and have all the general properties of 
spin-i fermion operators, even though additional quan- 
tum numbers (such as spin) are omitted since we consider 
only one band. 



B. Pairing term 

If the interband interaction is neglected, the interac- 
tion term is 



Tic 



i J2 I^2(ki,k2)ct_k,4ck,c_k., (21) 



ki,k2 



where T/2(ki,k2) = (*_ki*ki|V^2|*-k2*k2> is the two- 
particle matrix element. In the weak-coupling approachji 
one introduces the mean field potential 



A(k)=-^T/2(k,k')(ck'C_k' 



(22) 



and the interaction (|21|l is approximated by 

H2 = l ^[A(k)cl,cLk + A*(k)c_kCk], (23) 



where a constant term is neglected. It follows immedi- 
ately from the anticommutation of the fermion operators 
that 



A(-k) = -A(k) 



(24) 



One can derive the other transformation properties 
of A(k) from the fact that the full Hamiltonian H = 
Ti-i +Ti.2 is invariant under space group operations, time- 
reversal and gauge transformations. Using (|19|l . one ob- 
tains 

g : A(k) = A{g-'k) exp[-m<,(k) - m<,(-k)] (25) 

It will be seen already from a simple example given in 
Sec. IIVI that the phase factor in Eq. (|25|) is not triv- 
ial. In general, its dependence on g and k cannot be 
eliminated. Therefore, the function A(k) does not trans- 
form according to irreducible representations of the space 
(point) group. For crystals with a centre of symmetry, a 
comprehensive discussion of this property has been given 
by Blount i^ Instead of the pairing potential A^^ which 
transforms like {ipfj., ''py), an auxiliary object which trans- 
forms like {IC : tp^, ipy) is introduced. Here /i and v denote 
the set of single-particle quantum numbers. Following 
this idea, we define function x(k) by 



A(k) = x(k)i(k). 



(26) 



Then, using the definition of cil'^'' from the third Eq. 119|) 
the first term in H2 is expressed as 



i5:x(k)t(k)cJ,ct_, = l5]x(k)4cL( 



tec) 

k 



(27) 



Using the commutation property gtC = )Cg, it is straight- 



forward to show that 



10 



9 



JCC) 



exp[— iQ;g(k)]c 



t(^) 



Thus, from the invariance of 127|l . we find 
9 ■ X(k) = x(.9"^k). 



-Ik- 



(28) 



Eq. H28|) is the basis of the group-theoretical classifica- 
tion of SC states in noncentrosymmetric crystals. x(k) 
can be expanded in terms of basis functions Xi(k) of ir- 
reducible representations of the space (point) group. At 
this point we restrict our consideration to homogeneous 
SC states for simplicity, so that only the point group is 
involved. 



X(k) =^77jX*(k) 



(29) 



where iji will be identified as the components of the order 
parameter. It follows from Eqs. ifTTI) . (|^ and (ESI that 
x(— k) = x(k). Examples of even basis functions Xi(k) 
for the irreducible representations of point groups C^^ for 
CePtsSi and C2 for Ulr are given in Table Q] 

The other transformation properties of x(k) also follow 
from (O and ifTTj) . 

/C:x(k)=x*(k), C/i(0):x(k)=e2'^x(k), (30) 



TABLE I: Even basis functions for irreducible representations 
(IR) of point groups d^ and C2. The notations for IR and 
their character tables can be found in Ref. l2a . a, b, and c 
denote lattice constants. 



IR 


Non-periodic 


Periodic 


C4v 






Ai 


-L , K^ ~\~ f^iii •'j z 


cos kxa + cos kya, cos kzC 


A2 


f^xr^yyi^x "'!// 


sin fca;a sin kya{cos kxa — cos kyu) 


Bi 


i-2 i-2 


cos kxO. — cos kya 


B2 


i^xf^y 


sin kxasvakya 


E 


{ kxk. 


\ sin kx a sin kzC 




\ kyk. 


y smkyasinkzC 


C2 






^1 


r . f^x 1 '^v 1 '^ z 1 


cos kxa, cos kyb, cos kzC 




fCxIiy 


sin kxa sin kyb 


A2 


i^X^Z: f^yf^z 


sin kxa sin kzC, sin kyb sin kzC 



where we again use antilinearity of /C. As is usually done 
in the Ginzburg-Landau approach, the transformation 
properties of x(k) are reformulated as transformations 
of i]iJ If Xi(k) are chosen real, then g : rji = DikJik, 
^ ■ Vi = Vi s-iid Ui{9) : r]i = e'^^^rji, where D is the 
matrix corresponding to g in the representation chosen. 



C. Possible forms of the Ginzburg-Landau 
potential for CePtaSi 

Even though dv lacks inversion symmetry, the homo- 
geneous terms in the Ginzburg-Landau potential (GLP) 
coincide with those of D^hj^ because the GLP has to 
be invariant with respect to gauge transformations, in 
particular Ui{'k/2) : x(k) = — x(k). For every represen- 
tation F listed in Table Q] the product F (g) F* contains 
the representation Ai. The z-component of the gradient 
operator 



D 



^ 2ie . 

V A 

c 



(31) 



also transforms according to Ai, therefore the GLP den- 
sity contains a term linear in Z?^, 



For the two-dimensional representation E, 

F2D = Fo + J dr{Ki{\D,7^i\^ + \Dy7^2\^) 
+K2{\D,V2\^ + \DyTji\^) + K3[{D,r,i){Dyr,2)* + c.c] 
+K4[{D,T]2){Dym)* + c.c] -I- X5(|I?.ml' + \DzV2\'') 

2\2 



o*^2 



+a(lmr + \mn + /3i(hir + \m\'r + Mvlm - mri*2) 
+P,\rii?\m?}- 

(34) 



D. Green's functions 

The most straightforward way to calculate the normal 
G'(k, r) and anomalous i^(k, r), i^^(k, r) Green's func- 
tions within the weak-coupling approach is to use the 
Gor'kov equationsi22i2^ The Green's functions are defined 
as follows. 



G(k,r) = -(r.{ck(r)4(0)}) 

F(k,r) = (r,{ck(r)c_k(0)}) 

Ft(k,r) = (T.{cLk(r)4(0)}), 



(35) 



where T,- is the imaginary time ordering operator. The 
Gor'kov equations for the Fourier transforms are 






-ek]G(k,c^„) + A(k)^t(k^^^) ^ 1^ 
Ck]i^t(k,^„) + A*(k)G(k,L^„) = 0, (36) 



where a;„ — TTT{2n + 1) are the Matsubara frequencies 
for fermions. The equations are easily solved by 



G(k,c.„) = -(zc.„ + Ck)/(c^f,+ek + lx(k)n 
Ft(k,c.„) = x*(k)r(k)/(c^^^e^ + |x(k)n. 



(37) 



where we have made use of the fact that |A(k)p — 
lx(k)p. 

Gor'kov and Rashba^ considered a model of supercon- 
ductivity with split bands due to SOC and an isotropic 
pairing interaction. They showed that the thermodynam- 
ics of a such a superconductor is equivalent to that of a 
conventional s-wave superconductor. Here we have gen- 
eralized this result to the anisotropic case. As is evident 
from the denonrinator in Eq. I|37|l . the thernrodynamics 
of a noncentrosymnretric superconductor are governed by 
the gap |x(k)|. 



^VtDzTJi 



mDlvt- 



(32) 



III. JOSEPHSON EFFECT 



However, this terms is equal to the derivative -4- ^^ by^P- 
Hence, after integration over the sample, its contribution 
to the GLP vanishes. 

For all one-dimensional representations of Ci^,, the 
GLP is 

FiD = F^+ J dr[Kii\D,7^f + \Dy7^\^) + K2\D,7^f 

+a\r,\' + P\r^\% (33) 



We consider tunneling of SC electrons between a super- 
conductor with non-degenerate bands and a light conven- 
tional superconducting metal like Nb, in which electronic 
bands are spin degenerate. The tunneling Hamiltonian 
is 



Ti-T — 



E 

ki,k2,s 



Ts(ki,k2)c[^ak2,s+T;(ki,k2)a[^_^CkJ, 



(38) 



where, aj^ ^ creates an electron with spin s in the con- 
ventional superconductor and Ts(ki, k2) is the tunneling 
matrix element. By applying the time reversal operation 
to the first term in Ht, we obtain 



Y^ r;(ki,k2)i(ki)cLk,(-Jcry)5,/a_k2,s' 

ki,k2,s,s' 

= Yl ^(-ki)(Jo-j^)ss'r;,(-ki,-k2)4^ak, 

ki ,k2,s,s' 



(39) 



Therefore, 

T,(-ki,-k2) =t(ki)^(iay),,,T;,(ki,k2). (40) 



The following derivation follows that for conventional 
superconductors given in Refs. 13, l3.'i which is in turn 
based on the formalism originally proposed by Ambe- 
gaokar and BaratofF.'^'^ For simplicity, we restrict our con- 
sideration to the case of zero voltage. The current flowing 
between the two superconductors is j — e{dN{t)/dt) = 
ei{[nT,N]), where N{t) = e^^^* X^k CfcCke"*^^* is the 
time-dependent particle number operator. 



Following the same steps as described in Ref. | 
obtain 



= 2eTIm ^ T,(ki,k2)T,.(-ki, -k2)Ft(ki,c^„)^£"/(k2,cj„), 



(41) 



ki,k2, s,s',n 



where FfJ{'k,ujn) = t/'(k)(i(j)ss'/('^n + ^k + IV'(k)P) is the anomalous Green's function for a conventional supercon- 
ductor with the single-particle spectrum ek and singlet gap function ^(k).-^ Using (|37(l and H40|l . we obtain 



J = -2eTlm J^ E 



X*(ki)VXk2) 



7;i:, ^ K + e^, + |x(ki)P)(a.2 + el + |V'(k2)P) 

I 



El^^C'i'k^)!' 



(42) 



The sum over Matsubara frequencies can be easily calcu- 
lated if needed by using the formula 2T ^ .^ {u;f^+E^ )~^ — 
ta.nh{E / 2T) / E . Therefore, the final expression for the 
Josephson current does not depend on i(k). Once again, 
we see that the superconductor with SOC split bands 
behaves as a singlet superconductor with gap function 

X(k). 

The tunneling Hamiltonian Ht is invariant with re- 
spect to the operations g of the point group which leave 
the surface invariant. Thus, 



g : rs(ki,k2; 
It follows that 



— g-Mg(k 



'^V{g)ss'n'ig-'ku9-'k2 



(43) 



9- 



|T.(ki,k2) 



J2\Tsig-'k^,g-'k,)\' (44) 



Therefore, j is invariant under the operations g as ex- 
pected. 



IV. LIMIT OF SMALL SPIN-ORBIT COUPLING 

In order to verify the general results given in Sec. II, 
it is important to show that they remain valid for crys- 
tals with small SOC, i.e. for almost degenerate elec- 
tronic bands. In this Section we consider the Rashba 
approximation)^ because (i) the one particle Hamilto- 
nian is then exactly solvable and (ii) we can easily com- 
pare our results for superconductivity with those known 



from the literature, where the Rashba approximation was 
used.^4^ 

In the Rashba approximation, one assumes that the so- 
lutions of the Hamiltonian iJg = p^/2m + V{r) are plain 
waves with definite projections of spin and energy spec- 
trum e^, and that WV{r) can be replaced by a constant 
vector directed along the axis of symmetry n = (001). 
Therefore Eq. ^ is approximated by 



H[ = 



'P{ky + ikj;) 



jjyfvy IKx) 



(45) 



where /3 is a constant which characterizes the strength 
of the SOC. The Hamiltonian is diagonalized by the 



spmori 



.34 



* 



k.± 



iWr 



l^y \ ±iexp(i(/3k) 



(46) 



l|. Here, 
fc^. In the 



corresponding to the eigenvalues e^ = ^k ± /3|k 
exp(z(/Jk) = {k^^iky)l\\f.\\ and |k^| = Jk^. 
following, we shall consider coupling in '+' band as an 
example. 

In order to illustrate the complicated transformation 
properties H25|l of A(k), let us consider as an example 
g = m^;, the mirror plane perpendicular to x-axis. Acting 



by 771a; on ^k+, we obtain 




,z(k,m^r) 



,i(mj; k:,r) 



e^'^''*™,k,^47) 



where exp{iipm^k) = — exp(— i(^k) was used. Hence, 
the phase factor defined by Eq. 125|l is equal to {ikx + 

By applying time reversal to ^k,+ using Eq. (|f4|l . we 
find 

i(k) = iexp(-i</7k)- (48) 

The fermion operators corresponding to initial spin-up 
and spin-down states can be expressed in terms of the 
new band operators as 

1 

71 



"^kT - /^(^k+ + '^k-/ 



t ie-^ t t N 

c,'., = ^(Ck+-<-)- 



^ki 



V2 ^ ''+ 



(49) 



At this point we assume that before SOC has been turned 
on, the electrons were paired in the singlet state with gap 
function V'(k), 



7i^ = i^^(k)4^cLkx+h.c. 



(50) 



Using igni) and (gHl: we obtain 



^2 = ^E^wv'(k)4+c^-k+ + h. 



c. + 



(51) 



where the rest of the terms, which describe pairing in 
the other band and inter-band pairing, are neglected. By 
comparing (|51|l with (|23|) we obtain A(k) in the form 12()|) 
with x(k) = ^(k)/2. 

Interestingly, a similar result is obtained if one starts 
with triplet pairing for spin states^ 



^2 - 7;Y.^[~dx{^) + ^dy{\^)]ci^cl^^ + K(k) +zdj,(k)]4^cLk^ +4(k)(c[^clk4 + 4^cLkT)} +li.c. 



-Y,t{^)[dy{^)k - 4(k)fcy]c^cLk+ + li.c. 



(52) 



r 



where ki — fci/|k^| and again, only terms corresponding 
coupling within the '-!-' band are kept. Therefore, we 
again find A(k) in the form H26|l with 



H[: the singlet state transforms into itself as a scalar, 
whereas the triplet state transforms as a spin-1 statei^ 



X(k) = dy{\i)kx - d.j;(]i.)ky 



(53) 



Frigeri et al.^'' used weak coupling theory to show that 
a phase transition to the triplet pairing state d(k) = 
(— fcy, fca;, 0) (which transforms according to the repre- 
sentation A2) in CePtsSi is not suppressed by small 
but finite fi. According to l(5^. this corresponds to 
x(k) ex kl + ky, which is isotropic and fully gapped 
and transforms like Ai. By using the vectors d(k) cor- 
responding to the rest of the irreducible representations 
of Ciy (which are similar to those for D^h and are given 
in Ref. 7), the corresponding functions x(k) may be ob- 
tained from Eq. H53() and are listed in TablelTll Therefore, 
in general, starting from d(k) corresponding to a repre- 
sentation r, we generate x(k) using Eq. H53|l correspond- 
ing to another representation F', while in the singlet case, 
x(k) and V'(k)/2 belong to the same representation be- 
cause they are proportional. The difference is due to the 
transformation properties of the paired state under spin 
rotation, which was used to diagonalize the Hamiltonian 



TABLE II: Basis functions for triplet superconductivity in 
C4v crystals with non-degenerate bands and vanishing SOC 
and the corresponding even basis functions determined by 
Eq. ^. Note that A(k) is defined by Eq. ^ even for 
the triplet case. 



Al 


d(k) = (fc,,fc,,0) 


- 






d(k) = (fcS,fc^,0) 


X(k) = fc,fc^(fc^ - fc2) 


A2 


A2 


d(k) = (-fc^,fc,,0) 


x(k) = kl + kl 


Ai 


Si 


d(k) = (fc,,-fcj,,0) 


X(k) = k^ky 


B2 


B2 


d(k) = (fc^,fc,,0) 


x(k) = kl-kl 


Bi 


E 


di(k) = (fc,,0,0) 


Xi(k) = kyk^ 


E 




d2(k) = (0,fc,,0) 


X2(k) = k:,k^ 





DISCUSSION 



A close analogy between the superconductors under 
consideration and usual centrosymmetric singlet super- 
conductors should not be surprising. A particle described 
by the wave function (|10l) which is a solution of Hamil- 
tonian ^ has spin s(k) = ^(^'klcrj'J'k). Since \I>_k and 

^_J differ only by a constant phase factor, we have 



s(-k) = 



5< 
= /c 



'^".'1 



:(*kh 



rCC) 



*k) = 2(*k|<T 






l^k)' 



1 



= --(vl/kk|vl/k)=-s(k), 



(54) 



where we have used the fact that diagonal matrix el- 
ements of the spin operator are real and crC^^ — 
{—icFy)cr* [—iay]^^ — —a. Therefore, particles with op- 
posite momenta within the same band always have op- 
posite spin [See also Ref. [ij . 

It was first noted by Anderson^ that for crystals in 
which either the momentum k or spin is not a good quan- 
tum number, as in dirty superconductors, a one particle 
state should be paired with its time reversal. Using this 
idea, Samokhin tt alw^ define c_-^ as time reversal of cj^.. 
This in turn leads to the conclusion that A(k) is an odd 
function which transforms according to the irreducible 
representations of the point group. Then, the gaps in 
the quasiparticle spectrum of the superconducting states 
described by all one dimensional representations of C4U 
have line nodes. This, in particular, rules out the possi- 
bility of an isotropic full gap, and t herefore contradicts 
the theoretical results of Refs. I34l37l and the experimen- 
tal results for Cd2Re207. The origin of the discrepancy 
is that Anderson's statement is correct up to a phase fac- 
tor t(k), which is not important in the normal state due 
to gauge invariance but cannot be ignored in the super- 
conducting state42i In other words, one cannot define cj^ 
in an asymmetric unit of the Brillouin zone and then use 
the symmetry elements (including time reversal) to de- 
fine the states in the rest of the Brillouin zone. Such a 
procedure would lead to a discontinuity of the phase of 
the wave function on the boundaries of the asymmetric 
units. Also, if cj^. is paired with /CcJ^., then /CcJ^. must be 
paired with K?c\^ = — cj^, which is a contradiction. In- 
stead, if cj^ describes the particle with momentum k and 
k' = g~^k, where g is a point group element or time 
reversal, then the particle with momentum k' must be 
described by c! , — "' 



c^_i[^, which is proportional but not 



equal to g : c]^. 

The formalism developed in this paper can be di- 
rectly applied to Cd2Re207, which shows no sign of 
magnetic ordering and therefore possesses time reversal 
symmetry.^^ The low temperature structure has symme- 
try /4i22,^^ with point group D4, which is isomorphic to 
C41,. The superconducting order parameter corresponds 



to the representation Ai since no nodes in the quasi- 
particle spectrum have been found. 

In the application to CePtsSi, this theory should be 
slightly modified to include the effect of antiferromag- 
netic ordering (T/v = 2 K)iSi A neutron scattering study^ 
revealed that the antiferromagnetic structure is charac- 
terized by the wave vector Q = (0, 0, tt/c). Therefore, the 
normal state just above Tc is invariant under the opera- 
tion r(c)/C instead of /C, where t(c) is a lattice transla- 
tion along z-axis. This leaves valid all of the results of our 
discussion, since we consider translationally invariant su- 
perconducting states, and for the antiferromagnetic case, 
one still has e^ = £-k- However, the presence of antifer- 
romagnetic order can provide certain clues about possible 
superconducting state. It has been shown for zero SOC 
that the singlet superconducting state with gap function 
-(/•(k) is energetically favoured if ■(/'(k -I- Q) = — '0(k) [See 
Ref. and references therein]. Hence, the periodic basis 
functions in Table I, which depend on kz, may be good 
candidates. 

Even more exciting is the situation realized in UIr>S 
At ambient pressure it is ferromagnetic with Curie tem- 
perature of 46 K. Superconductivity occurs at high pres- 
sures of 2.6-2.7 GPa. One possibility is that the phase 
transition lines from the normal paramagnetic to ferro- 
magnetic state and from the normal paramagnetic to su- 
perconducting state meet at a quantum critical point. 
In that case superconducting states can again be clas- 
sified with respect to the symmetry of time reversal in- 
variant paramagnetic state. On the other hand, fcrro- 
magnetism can survive high pressures and coexist with 
the superconducting state. Then, time reversal symme- 
try is completely broken. Since the centre of symmetry is 
also missing, in general e\^ 7^ £-k- Therefore, a zero-field 
analog of the inhomogeneous Larkin-Ovchinnikov-Fulde- 
Ferrel state)^ may be realized in this material. On the 
other hand, superconductivity exists in a relatively nar- 
row pressure range. Hence, an alternative scenario may 
imply that e\^ and e_k are accidentally degenerate for a 
part of the Fermi surface. In either case, superconduc- 
tivity in Ulr deserves further investigation. 

In conclusion, we have considered the symmetry prop- 
erties of the gap function in superconductors with lifted 
spin degeneracy. We have shown that phase factors which 
appear with the gap function A(k) under point group op- 
erations and time reversal may be handled by the intro- 
duction of an even auxiliary function x(k), which trans- 
forms according to the irreducible representations of the 
point group. This function defines the behaviour of the 
superconductor in terms of thermodynamic and tunnel- 
ing properties. 

Note added. The criticism of Ref. [13 presented in this 
Article is addressed in Ref. [ij. 
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